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1 Introduction 

The synthesis of molecular magnets has undergone rapid progress in recent years HJEJEIEJEI 
|6). Each of the identical molecular units can contain as few as two and up to several dozens 
of paramagnetic ions ("spins"). One of the largest paramagnetic molecules synthesized to date, 
the polyoxometalate {Mo 72 Fe3o} contains 30 iron ions of spin s = 5/2. Although these 
materials appear as macroscopic samples, i. e. crystals or powders, the intermolecular mag- 
netic interactions are utterly negligible as compared to the intramolecular interactions. There- 
fore, measurements of their magnetic properties reflect mainly ensemble properties of single 
molecules. 

Their magnetic features promise a variety of applications in physics, magneto-chemistry, bi- 
ology, biomedicine and material sciences JH El El as well as in quantum computing ll8ll9lfT0l. 
The most promising progress so far is being made in the field of spin crossover substances using 
effects like "Light Induced Excited Spin State Trapping (LffiSST)" ifTTl . 
It appears that in the majority of these molecules the localized single-particle magnetic moments 
couple antiferromagnetically and the spectrum is rather well described by the Heisenberg model 
with isotropic nearest neighbor interaction sometimes augmented by anisotropy terms lfT2lfT3l 
[2J[T5J[I1|. Thus, the interest in the Heisenberg model, which is known already for a long time 
|[T71 . but used mostly for infinite one-, two-, and three-dimensional systems, was renewed by the 
successful synthesis of magnetic molecules. Studying such spin arrays focuses on qualitatively 
new physics caused by the finite size of the system. 

Several problems can be solved with classical spin dynamics, which turns out to provide accu- 
rate quantitative results for static properties, such as magnetic susceptibility, down to thermal 
energies of the order of the exchange coupling. However, classical spin dynamics will not be 
the subject of this chapter, it is covered in many publications on Monte-Carlo and thermostated 
spin dynamics. One overview article which discusses classical spin models in the context of 
spin glasses is given by Ref. lfT%1 . 

Theoretical inorganic chemistry itself provides several methods to understand and describe 
molecular magnetism, see for instance Ref. fl9]. In this chapter we would like to focus on 
those subjects which are of general interest in the context of this book. 

2 Substances 

From the viewpoint of theoretical magnetism it is not so important which chemical structures 
magnetic molecules actually have. Nevertheless, it is very interesting to note that they appear 
in almost all branches of chemistry. There are inorganic magnetic molecules like polyoxometa- 
lates, metal-organic molecules, and purely organic magnetic molecules in which radicals carry 
the magnetic moments. It is also fascinating that such molecules can be synthesized in a huge 
variety of structures extending from rather unsymmetric structures to highly symmetric rings. 
One of the first magnetic molecules to be synthesized was Mn-12-acetate ll20l (Mni 2 ) - 
[Mni20i2(CH 3 COO)i6(H 2 0)4] - which by now serves as the "drosophila" of molecular mag- 
netism, see e. g. dlZHSimiZH- A s shown in Fig.[l]the molecules contains four Mn(IV) ions 
(s = 3/2) and eight Mn(III) ions (s = 2) which are magnetically coupled to give an S = 10 
ground state. The molecules possesses a magnetic anisotropy, which determines the observed 
relaxation of the magnetization and quantum tunneling at low temperatures ll2Tll24l . 
Although the investigation of magnetic molecules in general - and of Mn-12-acetate in partic- 
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Fig. 1: Structure of Mn-12-acetate: On the l.h.s. the Mn ions are depicted by large spheres, on 
the r.h.s. the dominant couplings are given. With friendly permission by G. Chaboussant. 

ular - has made great advances over the past two decades, it is still a challenge to deduce the 
underlying microscopic Hamiltonian, even if the Hamiltonian is of Heisenberg type. Mn-12- 
acetate is known for about 20 years now and investigated like no other magnetic molecule, but 
only recently its model parameters could be estimated with satisfying accuracy ll25ll26l . 




Fig. 2: Structure of a chromium-4 and a chromium-8 ring. The Cr ions are depicted by large 
spheres. 

Another very well investigated class of molecules is given by spin rings among which iron rings 
("ferric wheels") are most popular fl27l EHl EH [30l EB [32l E21 E41. Iron-6 rings for instance 
can host alkali ions such as lithium or sodium which allows to modify the parameters of the 
spin Hamiltonian within some range ifTol I33L Another realization of rings is possible using 
chromium ions as paramagnetic centers. Figure |2] shows the structure of two rings, one with 
four chromium ions the other one with eight chromium ions. 

A new route to molecular magnetism is based on so-called Keplerate structures which allow the 
synthesis of truly giant and highly symmetric spin arrays. The molecule {Mo72Fe3o} 1171 1361 
containing 30 iron ions of spin s = 5/2 may be regarded as the archetype of such structures. 
Figure |3] shows on the l.h.s. the inner skeleton of this molecule - Fe and O-Mo-0 bridges - as 
well as the classical ground state [37 1 depicted by arrows on the r.h.s. OBI . 
One of the obvious advantages of magnetic molecules is that the magnetic centers of different 
molecules are well separated by the ligands of the molecules. Therefore, the intermolecular 
interactions are utterly negligible and magnetic molecules can be considered as being indepen- 
dent. Nevertheless, it is desirable to build up nano structured materials consisting of magnetic 
molecules in a controlled way. FigurelUgives an example of a planar structure consisting of lay- 
ers of {Mo7 2 Fe 3 o} ll3~Sll3l?ll which has been synthesized recently together with a linear structure 
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2.5 nm 




Fig. 3: Structure of {Mo^Fe^o}, a giant Keplerate molecule where 30 iron ions are placed at 
the vertices of an icosidodecahedron. L.h.s.: sketch of the chemical structure, r.h.s. magnetic 
structure showing the iron ions (spheres), the nearest neighbor interactions (edges) as well as 
the spin directions in the classical ground state. The dashed triangle on the l.h.s. corresponds 
to the respective triangle on the r.h.s.. With friendly permission by Paul Kdgerler H36V . 



Fig. 4: Square lattice of {Mo^Fe^} -molecules: Each molecule is connected with its four 
nearest neighbors by an antiferromagnetic coupling. With friendly permission by Paul Kdgerler 



consisting of chains of {Mo 72 Fe 3 o} ll40ll . These systems show new combinations of physical 
properties that stem from both molecular and bulk effects. 

Many more structures than those sketched above can be synthesized nowadays and with the 
increasing success of coordination chemistry more are yet to come. The final hope of course is 
that magnetic structures can be designed according to the desired magnetic properties. But this 
goal is not close at all, it requires further understanding of the interplay of magneto-chemistry 
and magnetic phenomena. One of the tools used to clarify such questions is density functional 
theory or other ab initio methods II3T1 1321 1331 1331 1351 1361. 

3 Theoretical techniques and results 

3.1 Hamiltonian 

It appears that in the majority of these molecules the interaction between the localized single- 
particle magnetic moments can be rather well described by the Heisenberg model with isotropic 
(nearest neighbor) interaction and an additional anisotropy term fTTl 031 IT31 IT51 IT6T Dipolar 
interactions are usually of minor importance. It is also found that antiferromagnetic interactions 
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are favored in most molecules leading to nontrivial ground states. 
Heisenberg Hamiltonian 

For many magnetic molecules the total Hamilton operator can be written as 

H -^Heisenberg anisotropy -^Zeeman (1) 

-^Heisenberg = ~ } u J U v£(u) • s(v) (2) 
u,v 
N 

^anisotropy = ~ ^ d u (e(u) ■ s(u)) 2 (3) 
«=1 

-^Zeeman = <?/^B-B • 5 . (4) 

The Heisenberg Hamilton operator 1 in the form given in Eq. © is isotropic, J uv is a symmet- 
ric matrix containing the exchange parameters between spins at sites u and v. The exchange 
parameters are usually given in units of energy, and J uv < corresponds to antiferromagnetic, 
J uv > to ferromagnetic coupling 2 . The sum in © runs over all possible tuples (u,v). The 
vector operators s(u) are the single-particle spin operators. 

The anisotropy terms © usually simplify to a large extend, for instance for spin rings, where 
the site-dependent directions e(u) are all equal, e. g. e(u) = e z and the strength as well is the 
same for all sites d u = d. 

The third part (Zeeman term) in the full Hamiltonian describes the interaction with the exter- 
nal magnetic field. Without singe-site and g-value anisotropy the direction of the field can be 
assumed to be along the z-axis which simplifies the Hamiltonian very much. 
Although the Hamiltonian looks rather simple, the eigenvalue problem is very often not solvable 
due to the huge dimension of the Hilbert space or because the number of exchange constants is 
too big to allow an accurate determination from experimental data. Therefore, one falls back 
to effective single- spin Hamiltonians for molecules with non-zero ground state spin and a large 
enough gap to higher-lying multiplets. 



Single-spin Hamiltonian 

For molecules like Mn 12 and Fe 8 which possess a high ground state spin and well separated 
higher lying levels the following single-spin Hamiltonian 

H = -D 2 S 2 Z - D A S A Z + H' (5) 
H' = gfi B B x S x (6) 

is appropriate, see e. g. Ref. ll23ll . The first two terms of the Hamilton operator H represent 
the anisotropy whereas H' is the Zeeman term for a magnetic field along the x-axis. The 
total spin is fixed, i. e. S — 10 for Mn 12 and Fe 8 , thus the dimension of the Hilbert space is 
dim(ft) =25 + 1. 

The effective Hamiltonian © is sufficient to describe the low-lying spectrum and phenomena 
like magnetization tunneling. Since H' does not commute with the ^-component of the total 

1 Operators are denoted by a tilde. 

2 One has to be careful with this definition since it varies from author to author 
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spin S z , every eigenstate | M } of S z , i. e. the states with good magnetic quantum number M, 
is not stationary but will tunnel through the barrier and after half the period be transformed into 

\-M). 



3.2 Evaluating the spectrum 

The ultimate goal is to evaluate the complete eigenvalue spectrum of the full Hamilton operator 
© as well as all eigenvectors. Since the total dimension of the Hilbert space is usually very 
large, e. g. dim(7i) = (2s + 1)^ for a system of N spins of equal spin quantum number s, a 
straightforward diagonalization of the full Hamilton matrix is not feasible. Nevertheless, very 
often the Hamilton matrix can be decomposed into a block structure because of spin symme- 
tries or space symmetries. Accordingly the Hilbert space can be decomposed into mutually 
orthogonal subspaces. Then for a practical evaluation only the size of the largest matrix to be 
diagonalized is of importance (relevant dimension). 



Product basis 

The starting point for any diagonalization is the product basis | m ) = | m 1; . . . , m ui . . . ,m N ) 
of the single-particle eigenstates of all s z (u) 

Sg(u) | mi, . . • , m u , . . . , m N ) = m u \ m u . . . , m u , . . . , m N ) . (7) 

These states are sometimes called Ising states. They span the full Hilbert space and are used to 
construct symmetry -related basis states. 



Symmetries of the problem 

Since the isotropic Heisenberg Hamiltonian includes only a scalar product between spins, this 
operator is rotationally invariant in spin space, i. e. it commutes with S and thus also with S z 



TT q2 

-£3 Heisenberg j & 



-^Heisenberg; $z 



0. (8) 



In a case where anisotropy is negligible a well-adapted basis is thus given by the simultaneous 
eigenstates | S, M, a ) of S 2 and S z , where a enumerates those states belonging to the same S 
and M HTIHSI . Since the applied magnetic field can be assumed to point into 2-direction for 
vanishing anisotropy the Zeeman term automatically also commutes with -^Heisenberg , S 2 , and 
S z . Since M is a good quantum number the Zeeman term does not need to be included in the 
diagonalization but can be added later. 

Besides spin symmetries many molecules possess spatial symmetries. One example is given by 
spin rings which have a translational symmetry. In general the symmetries depend on the point 
group of the molecule; for the evaluation of the eigenvalue spectrum its irreducible representa- 
tions have to be used [H31 IT6l HTll. Thus, in a case with anisotropy one looses spin rotational 
symmetries but one can still use space symmetries. Without anisotropy one even gains a further 
reduction of the relevant dimension. 
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Dimension of the problem 

The following section illuminates the relevant dimensions assuming certain symmetries 3 . 
If no symmetry is present the total dimension is just 

A' 



dim(W) = JJ (2s(«) + l) (9) 



u=l 



for a spin array of N spins with various spin quantum numbers. In many cases the spin quantum 
numbers are equal resulting in a dimension of the total Hilbert space of dim(7i) = (2s + 1)^. 
If the Hamiltonian commutes with S z then M is a good quantum number and the Hilbert space 
Tt can be divided into mutually orthogonal subspaces Tt(M) 



-Smax N 



H= H(M) , S max = s ( u ) • ( 10 ) 



M=-5 m ax u=l 



For given values of M, N and of all s(u) the dimension dim (H(M)) can be determined as 
the number of product states Q, which constitute a basis in H(M), with ^2 u m u = M. The 
solution of this combinatorial problem can be given in closed form [48 1 



dim (H(M)) 



(5 max -M)! 



^ \ Smax-M N ^ _ ^2s(»+l 



n 1 ^ 



dz . . 

3=1 J Z=0 



(ID 



For equal single-spin quantum numbers s(l) = • • • = s(N) = s, and thus a maximum total 
spin quantum number of S max = Ns, (fTTT) simplifies to 

d\m(H(M)) = f(N, 2s + 1, S max — M) with (12) 

m».)- E(- 1 )"(I)( Ar -rr" /1 )- 

n=0 \ / V / 

In both formulae (fTTTl and (TT2l . M may be replaced by \M\ since the dimension of TC(M) 
equals those of T~L(— M). \y / n\ in the sum symbolizes the greatest integer less or equal to z///i. 
Eq. (fT2l is known as a result of de Moivre ll49l . 

If the Hamiltonian commutes with S 2 and all individual spins are identical the dimensions of 
the orthogonal eigenspaces TC(S, M) can also be determined. The simultaneous eigenspaces 
H.(S, M) of S 2 and S z are spanned by eigenvectors of H. The one-dimensional subspace 
H(M = S max ) = H{S msai S^ax), especially, is spanned by | ), a state called magnon vacuum. 
The total ladder operators (spin rising and lowering operators) are 

S ± = S x ±iS y . (13) 

For S > M, S~ maps any normalized iJ-eigenstate G H(S, M + 1) onto an i/-eigenstate 
G H(S, M) with norm a/ S(S + 1) — M(M + 1). 

For < M < TC(M) can be decomposed into orthogonal subspaces 

H(M) = H(M, M) © S~H(M + 1) (14) 



Work done with Klaus Barwinkel and Heinz-Jiirgen Schmidt, Universitat Osnabriick, Germany. 
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with 

S~H(M+1) = H(S,M). (15) 

S>M+1 

In consequence, the diagonalization of H in H has now been traced back to diagonalization in 
the subspaces H(S, S), the dimension of which are for S < S m:a 

dim (H(S, S)) = dim (H(M = S)) - dim (H(M = S + 1)) (16) 

and can be calculated according to (fT2l . 

As an example for space symmetries I would like to discuss the translational symmetry found 
in spin rings. The discussed formalism can as well be applied to other symmetry operations 
which can be mapped onto a translation. Any such translation is represented by the cyclic shift 
operator T or a multiple repetition. T is defined by its action on the product basis © 

T | mi, ... , mjv-i, m N ) = \ m N , mi, ... , mjy-i ) • (IV) 
The eigenvalues of T are the iV-th roots of unity 

z fc = exp j-z^J , fc = 0,...,JV-l, p k = 27rk/N, (18) 

where k will be called translational (or shift) quantum number and p k momentum quantum 
number or crystal momentum. The shift operator T commutes not only with the Hamiltonian but 
also with total spin. Any H(S, M) can therefore be decomposed into simultaneous eigenspaces 
H{S, M,k) of | 2 , S z and T. 

In the following we demonstrate how an eigenbasis of both S z and T can be constructed, this 
basis spans the orthogonal Hilbert spaces TC(M, k). How total spin can be included by means 
of an irreducible tensor operator approach is described in Refs. |[T31[T6lH7l . 
A special decomposition of H into orthogonal subspaces can be achieved by starting with the 
product basis and considering the equivalence relation 

10)^ |V)=ri0>» ne{l,2,...,iV} (19) 

for any pair of states belonging to the product basis. The equivalence relation then induces 
a complete decomposition of the basis into disjoint subsets, i. e. the equivalence classes. A 
"cycle" is defined as the linear span of such an equivalence class of basis vectors. The obviously 
orthogonal decomposition of TC into cycles is compatible with the decomposition of H into the 
various TC(M). Evidently, the dimension of a cycle can never exceed N. Cycles are called 
"proper cycles" if their dimension equals N, they are termed "epicycles" else. One of the iV 
primary basis states of a proper cycle may arbitrarily be denoted as 

\ifa)= \ m 1 ,...,m N ) (20) 

and the remaining ones may be enumerated as 

|Vwi>=r I -01 > , n = l,2,...,AT-l. (21) 
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The cycle under consideration is likewise spanned by the states 

l^) = -^E O^'llM (22) 

which are eigenstates of T with the respective shift quantum number k. Consequently, every 
k occurs once in a proper cycle. An epicycle of dimension D is spanned by D eigenstates of 
T with each of the translational quantum numbers k = 0, N/D, . . . , (D — 1)N/D occurring 
exactly once. 

As a rule of thumb one can say that the dimension of each H(M, k) is approximately dim(7Y(M, k)) « 
dirn(W(M))/iV. An exact evaluation of the relevant dimensions for spin rings can be obtained 
from Ref. B8I . 



Exact diagonalization 

If the relevant dimension is small enough the respective Hamilton matrices can be diagonalized, 
either analytically ll50l ISTl 1481 or numerically, see e. g. ll5^l5^l54ll55lfT3ll56ll57ll47l . 
Again, how such a project is carried out, will be explained with the help of an example, a simple 
spin ring with N = 6 and s = 1/2. The total dimension is dim(7i) = (2s + 1)^ = 64. The 
Hamilton operator © simplifies to 

N 

#Heisenberg = ~2J ^ s( U ) ■ s( U + 1) , N + 1 = 1 . (23) 

u=l 

We start with the magnon vacuum \VL) = | + + + + ++) which spans the Hilbert space 
H(M) with M = Ns = 3. "±" are shorthand notations for m = ±1/2. The dimension of the 
subspace dim(7Y(M = Ns)) is one and the energy eigenvalue is Eq = —2JNs 2 = —3 J. | Vl ) 
is an eigenstate of the shift operator with k = 0. Since S is also a good quantum number in this 
example | Q ) has to be an eigenstate of S 2 , too, the quantum number is S = Ns. 
The next subspace H(M) with M = Ns — 1 = 2 is spanned by | — h + + ++ ) and the five 
other vectors which are obtained by repetitive application of T. This subspace obviously has 
the dimension N, and the cycle spanned by T n \ — h + + ++ ), n = 0, . . . , N — 1 is a proper 
one. Therefore, each k quantum number arises once. The respective eigenstates of T can be 
constructed according to Eq. (l2"2l as 

\M = 2,k) = J2 (e^lY {- + + + ++) . (24) 

All subspaces Ti.(M, k) have dimension one. Since S~ \ £1 ) is a state belonging to H(M = 
Ns — 1) with the same fc-quantum number as [ O ) it is clear that | M = 2, k = ) is a already 
an eigenstate of S 2 with S = Ns. The other | M = 2, k ^ ) must have S = Ns - 1. 
The next subspace H(M) with M = Ns — 2 = 1 is spanned by three basic vectors, i. e. 

| h + ++), | — I h ++), | — h H h+) and the repetitive applications of T onto 

them. The first two result in proper cycles, the third vector | — h H h+ ) results in an epicy- 
cle of dimension three, thus for the epicycle we find only k quantum numbers k = 0,2, 4. The 
energy eigenvalues found in the subspace H(M = Ns — 1) ("above") must reappear here which 
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again allows to address an S quantum number to these eigenvalues. The dimension of the sub- 
space TC(M = 1) is 15, the dimensions of the subspaces H(M, k) are 3 (k — 0), 2 (k = 1), 3 
(k = 2), 2(k — 3), 3 (k = 4), and 2(k — 5). 

The last subspace which has to be considered belongs to M = and is spanned by | h ++ ) , | 1 h4 

and repetitive applications of T. Its dimension is 20. Here | — I 1 h ) leads to an epicycle 

of dimension two. 

The Hamilton matrices in subspaces with M < need not to be diagonalized due to the S z - 
symmetry, i. e. eigenstates with negative M can be obtained by transforming all individual 
m u — > — m u . Summing up the dimensions of all H(M) then yields 1 + 6 + 15 + 20 + 15 + 6 + 1 = 

64 y/. 





Fig. 5: Energy eigenvalues as a function of total spin quantum number S (l.h.s.) and k (r.h.s.). 

Figure |5] shows the resulting energy spectrum both as a function of total spin S as well as a 
function of translational quantum number k. 



Projection and Lanczos method 

Complex hermitian matrices can be completely diagonalized numerically up to a size of about 
10,000 by 10,000 which corresponds to about 1.5 Gigabyte of necessary RAM. Nevertheless, 
for larger systems one can still use numerical methods to evaluate low-lying energy levels and 
the respective eigenstates with high accuracy. 

A simple method is the projection method 11531 which rests on the multiple application of the 
Hamiltonian on some random trial state. 

To be more specific let's approximate the ground state of a spin system. We start with a random 
trial state | O ) an d apply an operator which "cools" the system. This operator is given by the 
time evolution operator with imaginary time steps 

\4>i) = exp{- £ #} \ <j> ). (25) 

Expanding | O ) mt0 eigenstates | v ) of the Hamilton operator elucidates how the method 
works 

= J2 ex v{~ £E »} IWH0o> (26) 
= exp{- £ £ }^exp -£ )} \v)(v\fa) . (27) 

Relative to the ground state component all other components are exponentially suppressed. For 



Quantum Theory of Molecular Magnetism 



X5.ll 



practical purposes equation (l26l) is linearized and recursively used 



= 11-eH) {fa) , |fr +1 ) = — j j <+1 > (28) 



( <l>i+l I <j>i+l ) 

e has to be small enough in order to allow the linearization of the exponential. It is no problem 
to evaluate several higher-lying states by demanding that they have to be orthogonal to the 
previous ones. Restricting the calculation to orthogonal eigenspaces yields low-lying states in 
these eigenspaces which allows to evaluate even more energy levels. The resulting states obey 
the properties of the Ritz variational principle, i. e. they lie above the ground state and below 
the highest one. 

Another method to partially diagonalize a huge matrix was proposed by Cornelius Lanczos in 
1950 ll58ll59l . Also this method uses a (random) initial vector. It then generates an orthonormal 
system in such a way that the representation of the operator of interest is tridiagonal. Every 
iteration produces a new tridiagonal matrix which is by one row and one column bigger than 
the previous one. With growing size of the matrix its eigenvalues converge against the true ones 
until, in the case of finite dimensional Hilbert spaces, the eigenvalues reach their true values. 
The key point is that the extremal eigenvalues converge rather quickly compared to the other 
ones [ 60 1 . Thus it might be that after 300 Lanczos steps the ground state energy is already 
approximated to 10 figures although the dimension of the underlying Hilbert space is 10 8 . 
A simple Lanczos algorithm looks like the following. One starts with an arbitrary vector | ip ), 
which has to have an overlap with the (unknown) ground state. The next orthogonal vector is 
constructed by application of H and projecting out the original vector | ip ) 

|Vi) = (l- \ 4> HA\)H\1> ) = H\4> )-(i!>o\H\il> )\iJ> ), (29) 
which yields the normalized vector 

|^)= ; =. (30) 

Similarly all further basis vectors are generated 

IVk+i) = (1- l^fcXVfcl - \ip k -i)(^k-x\)H\^ k ) (31) 
= H | ip k ) - ( ip k | H | ip k ) | ip k ) - ( | H | i[> k ) | ip k -i ) 

and 

\M= 7=^4 r (32) 

The new Lanczos vector is by construction orthogonal to the two previous ones. Without proof 
we repeat that it is then also orthogonal to all other previous Lanczos vectors. This constitutes 
the tridiagonal form of the resulting Hamilton matrix 

Tij = ( ipi | H | ipj ) with T itj = if \i - j\ > 1 . (33) 



The Lanczos matrix T can be diagonalized at any step. Usually one iterates the method until a 
certain convergence criterion is fulfilled. 
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The eigenvectors of H can be approximated using the eigenvectors | M } of T 

n 

(34) 

i=0 

where /i labels the desired energy eigenvalue, e. g. the ground state energy, n denotes the 
number of iterations. 

The simple Lanczos algorithm has some problems due to limited accuracy. One problem is that 
eigenvalues may collapse. Such problems can be solved with more refined formulations of the 
method ll59l . 

DMRG 

The DMRG technique lloTll has become one of the standard numerical methods for quantum 
lattice calculations in recent years ll62l 1631 . Its basic idea is the reduction of Hilbert space 
while focusing on the accuracy of a target state. For this purpose the system is divided into 
subunits - blocks - which are represented by reduced sets of basis states. The dimension m of 
the truncated block Hilbert space is a major input parameter of the method and to a large extent 
determines its accuracy. 

DMRG is best suited for chain-like structures. Many accurate results have been achieved by 
applying DMRG to various (quasi-)one-dimensional systems Il64ll56l f65l. The best results were 
found for the limit of infinite chains with open boundary conditions. It is commonly accepted 
that DMRG reaches maximum accuracy when it is applied to systems with a small number of 
interactions between the blocks, e. g. systems with only nearest-neighbor interaction ll62ll . 
It is not a priori clear how good results for finite systems like magnetic molecules are 4 . Such 
systems are usually not chain-like, so in order to carry out DMRG calculations a mapping onto 
a one-dimensional structure has to be performed [62|. Since the spin array consists of a count- 
able number of spins, any arbitrary numbering is already a mapping onto a one-dimensional 
structure. However, even if the original system had only nearest-neighbor exchange, the new 
one-dimensional system has many long-range interactions depending on the way the spins are 
enumerated. Therefore, a numbering which minimizes long range interactions is preferable. 
Fig.|6]shows the graph of interactions for the molecule {Mo 72 Fe3 } which we want to consider 
as an example in the following ll66ll . 




Fig. 6: One-dimensional projection of the icosidodecahedron: the lines represent interactions. 

For finite systems a block algorithm including sweeps, which is similar to the setup in White's 
original article lloTI . has turned out to be most efficient. Two blocks are connected via two 
single spin sites, these four parts form the superblock see Fig. IT] 

For illustrative purposes we use a simple Heisenberg Hamiltonian, compare ©. The Hamilto- 
nian is invariant under rotations in spin space. Therefore, the total magnetic quantum number 

4 Work done with Matthias Exler, Universitat Osnabruck, Germany. 
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Fig. 7: Block setup for DMRG "sweep" algorithm: The whole system of N spins constitutes 
the superblock. The spins {1,2, ... ,p} belong to the left block, the other spins {p+1, ... , N} 
to the right block. 




Fig. 8: Two-dimensional projection of the icosidodecahedron, the site numbers are those used 
in our DMRG algorithm. 

M is a good quantum number and we can perform our calculation in each orthogonal subspace 
H(M) separately. 

Since it is difficult to predict the accuracy of a DMRG calculation, it is applied to an exactly 
diagonalizable system first. The most realistic test system for the use of DMRG for {Mo 72 Fe3o} 
is the icosidodecahedron with spins s = 1/2. This fictitious molecule, which possibly may be 
synthesized with vanadium ions instead of iron ions, has the same structure as {Mo7 2 Fe 30 }, 
but the smaller spin quantum number reduces the dimension of the Hilbert space significantly. 
Therefore a numerically exact determination of low-lying levels using a Lanczos method is 
possible [67|. These results are used to analyze the principle feasibility and the accuracy of the 
method. 

The DMRG calculations were implemented using the enumeration of the spin sites as shown in 
Figs. |^1 and [S] This enumeration minimizes the average interaction length between two sites. 
In Fig. 13 the DMRG results (crosses) are compared to the energy eigenvalues (circles) deter- 
mined numerically with a Lanczos method 115711551 . Very good agreement of both sequences, 
with a maximal relative error of less than 1% is found. Although the high accuracy of one- 
dimensional calculations (often with a relative error of the order of 10 -6 ) is not achieved, the 
result demonstrates that DMRG is applicable to finite 2D spin systems. Unfortunately, increas- 
ing m yields only a weak convergence of the relative error, which is defined relative to the width 
of the spectrum 

/ \ -&DMRG (m) - Eq 

€(m)= \E**-E*\ ■ (35) 

The dependence for a quasi two-dimensional structure like the icosidodecahedron is approx- 
imately proportional to 1/to (see Fig. flOb. Unfortunately, such weak convergence is charac- 
teristic for two-dimensional systems in contrast to one-dimensional chain structures, where the 
relative error of the approximate energy was reported to decay exponentially with m (6TJ. Nev- 
ertheless, the extrapolated ground state energy for s = 1/2 deviates only by e = 0.7 % from the 
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Fig. 9: Minimal energy eigenvalues of the s = 1/2 icosidodecahedron. The DMRG result with 
m = 60 is depicted by crosses, the Lanczos values by circles. The rotational band is discussed 
in subsection \3.4\ 

ground state energy determined with a Lanczos algorithm. 



□ 8=1/2 



0.01 0.015 
1/m 



Fig. 10: Dependence of the approximate ground state energy on the DMRG parameter m. E 
is the true ground state energy in the case s = 1/2 and the extrapolated one for s = 5/2. 



The major result of the presented investigation is that the DMRG approach delivers acceptable 
results for finite systems like magnetic molecules. Nevertheless, the accuracy known from one- 
dimensional systems is not reached. 



Spin-coherent states 

Spin-coherent states [ 68 ] provide another means to either treat a spin system exactly and inves- 
tigate for instance its dynamics [69| or to use spin coherent states in order to approximate the 
low-lying part of the spectrum. They are also used in connection with path integral methods. In 
the following the basic ideas and formulae will be presented. 

The obvious advantage of spin-coherent states is that they provide a bridge between classical 
spin dynamics and quantum spin dynamics. Spin coherent states are very intuitive since they 
parameterize a quantum state by the expectation value of the spin operator, e. g. by the two 
angles which represent the spin direction. 
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Spin coherent states | z ) are defined as 

\z) = jn ^ s f: ] l^)^\s,rn = s-p), zeC. (36) 

In this definition spin-coherent states are characterized by the spin length s and a complex value 
z. The states (I36T) are normalized but not orthogonal 

/ i \ 1 / i \ (l + y*z) 2s 

(z \ z) = 1, ly \ z) — -, , . , , n . . 

(1 + M 2 ) s (i + \A 2 Y 

Spin-coherent states provide a basis in single-spin Hilbert space, but they form an overcomplete 
set of states. Their completeness relation reads 



2S + 1 [ A*~ \ Z )( Z 



1 = — / d V ' ; , ' , d 2 z = dRe(z) dlm(z) . (37) 

~ 7T / (1 + \Z\ 2 Y 



The intuitive picture of spin-coherent states becomes obvious if one transforms the complex 
number z into angles on a Riemann sphere 

z = tan(0/2)e** , < 9 < it , O<0<2tt. (38) 

Thus, spin-coherent states may equally well be represented by two polar angles 9 and </>. Then 
the expectation value of the spin operator s is simply 

(sin(0) cos(0)\ 
sm(9) sin(0) . (39) 
cos(£) / 

Using d3*Sb the definition of the states | 9, <\> ) which is equivalent to Eq. (13*61 ) is then given by 

= y(^) [ cos (^/ 2 )] (2 "~ P) [e^sin(0/2)] p \s,m = s-p) (40) 
and the completeness relation simplifies to 

I = J dfi \9,<P)(0A\ ■ (41) 

Product states of spin-coherent states span the many-spin Hilbert space. A classical ground 
state can easily be translated into a many-body spin-coherent state. One may hope that this 
state together with other product states can provide a useful set of linearly independent states in 
order to approximate low-lying states of systems which are too big to handle otherwise. But it 
is too early to judge the quality of such approximations. 

3.3 Evaluation of thermodynamic observables 

For the sake of completeness we want to outline how basic observables can be evaluated both 
as function of temperature T and magnetic field B. We will assume that H, S z 



for this 
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part, so that the energy eigenvectors | v ) can be chosen as simultaneous eigenvectors of S z 
with eigenvalues E V (B) and M v . The energy dependence of E P (B) on B is simply given by 
the Zeeman term. If H and S z do not commute the respective traces for the partition function 
and thermodynamic means have to be evaluated starting from their general definitions. 
The partition function is defined as follows 

Z(T,B) = tr[ e - %} = J2 e ~^ iB) ■ ^ 

V 

Then the magnetization and the susceptibility per molecule can be evaluated from the first and 
the second moment of S z 

M(T,B) = -Itr{s/i B 5 z e-^} (43) 



X(T,B) = (44) 



k B T 



2 ^z2 M ^-" EAB) -(^zZ M »> 



-PE V (B) 



In a similar way the internal energy and the specific heat are evaluated from first and second 
moment of the Hamiltonian 

U(T,B) = -Atrjtfe-^} (45) 



-\Y, e ^ b ) q ~ peab) 



C(T,B) = E2±*L (46) 



Z 
OB 

k B 



■PE„(B) 



3.4 Properties of spectra 

In the following chapter I am discussing some properties of the spectra of magnetic molecules 
with isotropic and antiferromagnetic interaction. 



Non-bipartite spin rings 

With the advent of magnetic molecules it appears to be possible to synthesize spin rings with 
an odd number of spins. Although related to infinite spin rings and chains such systems have 
not been considered mainly since it does not really matter whether an infinite ring has an odd or 
an even number of spins. In addition the sign rule of Marshall and Peierls lITOll and the famous 
theorems of Lieb, Schultz, and Mattis ll7Tll72ll provided valuable tools for the understanding of 
even rings which have the property to be bipartite and are thus non-frustrated. These theorems 
explain the degeneracy of the ground states in subspaces TC(M) as well as their shift quantum 
number k or equivalently crystal momentum quantum number pk = 2nk/N. 
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3 


4 


5 


6 


7 


8 


9 


10 






1.5 


0.5 


1 


0.747 


0.934 


0.816 


0.913 


0.844 


0.903 


E /{NJ) 


1 

Z 


1 


4 


1 


4 


1 


4 


1 


4 


1 


deg 







1/2 





1/2 





1/2 





1/2 





S 




1 


1,2 





1,4 


3 


2,5 





2,7 


5 


k 




4.0 


3.0 


2.0 


2.236 


1.369 


2.098 


1.045 


1.722 


0.846 


AE/\J\ 
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o 
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3 
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3/2 
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1/2 
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3/2 
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3/2 
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S 










2 








1,6 
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3,6 





k 




4 


2 


3 


2.612 


2.872 


2.735 


2.834 


2.773 


2.819 


E /(NJ) 
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1 


1 
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S 































k 




4.0 


2.0 


2.0 


1.929 


1.441 


1.714 


1.187 


1.540 


1.050 


AE/\J\ 


1 


3 


9 


3 


6 


3 


6 


3 


6 


3 


deg 




1 


1 


1 


1 


1 


1 


1 


1 


1 


S 




1 


0,1,2 


2 


2,3 


3 


3,4 


4 


4,5 


5 


k 



Table 1: Properties of ground and first excited state of AF Heisenberg rings for various N and 
s: ground state energy E , gap AE, degeneracy deg, total spin S and shift quantum number k. 

Nowadays exact diagonalization methods allow to evaluate eigenvalues and eigenvectors of H 
for small even and odd spin rings of various numbers N of spin sites and spin quantum numbers 
s where the interaction is given by antiferromagnetic nearest neighbor exchange Il52ll53ll54ll73l 
|74l|75]|. Although Marshall-Peierls sign rule and the theorems of Lieb, Schultz, and Mattis do 
not apply to non-bipartite rings, i. e. frustrated rings with odd N, it turns out that such rings 
nevertheless show astonishing regularities 5 . Unifying the picture for even and odd N, we find 
for the ground state without exception 117411751 : 

1 . The ground state belongs to the subspace H(S) with the smallest possible total spin quan- 
tum number S; this is either S = for N-s integer, then the total magnetic quantum 
number M is also zero , or S = 1 / 2 f or N • s half integer, then M = ±1/2. 

2. If N-s is integer, then the ground state is non-degenerate. 

3. If N-s is half integer, then the ground state is fourfold degenerate. 

4. If s is integer or N-s even, then the shift quantum number is k = 0. 

5. If s is half integer and N-s odd, then the shift quantum number turns out to be k = N/2. 

6. If N-s is half integer, then k = [(N + 1)/4J and k = N — [(N + 1)/4J is found. 
L(iV + 1)/4J symbolizes the greatest integer less or equal to (N + l)/4. 

5 Work done with Klaus Barwinkel and Heinz-Jiirgen Schmidt, Universitat Osnabriick, Germany. 
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1.411 
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1.838 


0.938++ 


AE/\J\ 


3 
2 


3 


16 


3 


8 


3 


8 


3 


8 


3 


deg 




1 


3/2 


1 


3/2 


1 


3/2 


1 


3/2 


1 


S 







0, 1,2 


2 


2,3 





1,6 


4 


3,6 





k 




12 


6 


10 


8.456 


9.722 


9.045 


9.630 


9.263 tf 


9.590++ 


E /{NJ) 


2 


1 


1 


1 


1 


1 


1 


1 


1 


1 


deg 































S 































k 




4.0 


2.0 


2.0 


1.922 


1.394 


1.652 


1.091 


1.431*+ 


0.906*+ 


AE/\J\ 

I I I 


2 


3 


9 


3 


6 


3 


6 


3 


6 


3 


deg 




1 


1 


1 


1 


1 


1 


1 


1 


1 


S 




1 
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14.528 t 


13.848 tf 


14.475*+ 


E /(NJ) 


5 
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1 


4 
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1/2 





1/2 





1/2 
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1,2 
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Table 2: Properties of ground and first excited state of AF Heisenberg rings for various N 
and s (continuation): ground state energy E , gap AE, degeneracy deg, total spin S and shift 
quantum number k. f - 0. Waldmann, private communication, f f - projection method [55]. 

In the case of s = 1/2 one knows the /c-quantum numbers for all iV via the Bethe ansatz 
ll54ll73ll . and for spin s = 1 and even N the k quantum numbers are consistent with Ref. ll53ll . 
It appears that for the properties of the first excited state such rules do not hold in general, but 
only for "high enough" N > 5 11731. Then, as can be anticipated from tables [l] and El we can 
conjecture that 

• if iV is even, then the first excited state has S = 1 and is threefold degenerate, and 

• if N is odd and the single particle spin is half-integer, then the first excited state has 
S = 3/2 and is eightfold degenerate, whereas 

• if iV is odd and the single particle spin is integer, then the first excited state has S = 1 
and is sixfold degenerate. 

Considering relative ground states in subspaces H(M) one also finds - for even as well as for 
odd N - that the shift quantum numbers k show a strikingly simple regularity for N ^ 3 

N 

k = ±(Ns - M) \— 1 mod N , (47) 



Quantum Theory of Molecular Magnetism 



X5.19 



where \N/2~\ denotes the smallest integer greater than or equal to N/2 ll76ll . For N = 3 
and 3s — 2 > \M\ > 1 one finds besides the ordinary A>quantum numbers given by d47t 
extraordinary A>quantum numbers, which supplement the ordinary ones to the complete set 

{A;} = {0, 1, 2}. 

For even N the k values form an alternating sequence 0, N/2, 0, N/2, ... on descending from 
the magnon vacuum with M = Ns as known from the sign-rule of Marshall and Peierls lITOll . 
For odd N it happens that the ordinary A>numbers are repeated on descending from M < Ns—1 
to M - 1 iff N divides [2(Ns — M) + 1]. 

Using the A; -rule one can as well derive a rule for the relative ground state energies and for the 
respective S quantum numbers: 

• For the relative ground state energies one finds that if the A>number is different in adjacent 
subspaces, E min (S) < E min (S + 1) holds. If the fc-number is the same, the energies could 
as well be the same. 

• Therefore, if iV (even or odd) does not divide (2(Ns — M) + 1) \N/2~\ , then any relative 
ground state in TC(M) has the total spin quantum number S — \M\. 

• This is always true for the absolute ground state which therefore has S = for Ns integer 
and 5 = 1/2 for Ns half integer. 

The £>rule d47t is founded in a mathematically rigorous way for N even [17011711 [721 . N = 3, 
M = Ns, M = Ns — 1, and M = Ns — 2 iTToll . An asymptotic proof for large enough iV can 
be provided for systems with an asymptotically finite excitation gap, i. e. systems with integer 
spin s for which the Haldane conjecture applies 11771 1781. In all other cases numerical evidence 
was collected and the A; -rule as a conjecture still remains a challenge 1761 . 

Rotational bands 

For many spin systems with constant isotropic antiferromagnetic nearest neighbor Heisenberg 
exchange the minimal energies E min (S) form a rotational band, i. e. depend approximately 
quadratically on the total spin quantum number S 11791 18U1I8T1 

D(N s) 

E mm (S) ~E a -J S(S + 1) . (48) 

The occurrence of a rotational band has been noted on several occasions for an even number of 
spins defining a ring structure, e. g. see Ref. lUffl . The minimal energies have been described as 
"following the Lande interval rule" Il28ll29ll30ll32ll . However, we 6 find that the same property 
also occurs for rings with an odd number of spins as well as for the various polytope config- 
urations we have investigated, in particular for quantum spins positioned on the vertices of a 
tetrahedron, cube, octahedron, icosahedron, triangular prism, and an axially truncated icosa- 
hedron. Rotational modes have also been found in the context of finite square and triangular 
lattices of spin- 1/2 Heisenberg antiferromagnets ll82ll83l . 

There are several systems, like spin dimers, trimers, squares, tetrahedra, and octahedra which 
possess a strict rotational band since their Hamiltonian can be simplified by quadrature. As an 



6 Work done together with Marshall Luban, Ames Lab, Iowa, USA. 
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example the Heisenberg square, i. e., a ring with N = 4 is presented. Because the Hamilton 
operator (l23l) can be rewritten as 

H = —J (s 2 - S 2 13 - g 4 ) , (49) 
| 13 = |(1) + |(3) , £ 24 = |(2) + |(4) , (50) 

with all spin operators S 2 , S 2 3 and S 2 A commuting with each other and with H, one can directly 
obtain the complete set of eigenenergies, and these are characterized by the quantum numbers 
S, S\ 3 and S 2 a- In particular, the lowest energy for a given total spin quantum number S occurs 
for the choice £13 = S24 = 2s 

E min (S) = -J [5 (5 + 1) - 2 • 2s (2s + 1)] = Eq- J 5 (5 + 1) , (51) 

where E = 4s (2s + 1) J is the exact ground state energy. The various energies E min (S) form 
a rigorous parabolic rotational band of excitation energies. Therefore, these energies coincide 
with a parabolic fit (crosses connected by the dashed line on the l.h.s. of Fig. [Til passing 
through the antiferromagnetic ground state energy and the highest energy level, i. e., the ground 
state energy of the corresponding ferromagnetically coupled system. 




Fig. 11: Energy spectra of antiferromagnetically coupled Heisenberg spin rings (horizontal 
dashes ). The crosses connected by the dashed line represent the fit to the rotational band ac- 
cording to (I5TT) . which matches both the lowest and the highest energies exactly. On the l.h.s the 
dashed line reproduces the exact rotational band, whereas on the r.h.s. it only approximates it, 
but to high accuracy. The solid line on the r.h.s. corresponds to the approximation Eq. ( 152 D . 



It turns out that an accurate formula for the coefficient D(N, s) of (BTT) can be developed using 
the sublattice structure of the spin array J79). As an example we repeat the basic ideas for 
Heisenberg rings with an even number of spin sites ll32ll . Such rings are bipartite and can be 
decomposed into two sublattices, labeled A and B, with every second spin belonging to the 
same sublattice. The classical ground state (Neel state) is given by an alternating sequence 
of opposite spin directions. On each sublattice the spins are mutually parallel. Therefore, a 
quantum trial state, where the individual spins on each sublattice are coupled to their maximum 
values, Sa = Sb = Ns/2, could be expected to provide a reasonable approximation to the true 
ground state, especially if s assumes large values. For rings with even iV the approximation to 
the respective minimal energies for each value of the total spin S = Sa + Sb is then given by 
E3 



E 



approx 

min 



4 J 
'~N 



S{S + l)-2 



Ns [Ns 



+ 1 



(52) 



Quantum Theory of Molecular Magnetism 



X5.21 



This approximation exactly reproduces the energy of the highest energy eigenvalue, i. e., the 
ground state energy of the corresponding ferromagnetically coupled system (S = Ns). For all 
smaller S the approximate minimal energy E^™ X (S) is bounded from below by the true one 
(Rayleigh-Ritz variational principle). The solid curve displays this behavior for the example of 
iV = 6, s = 3/2 in Fig.[TT](r.h.s.). The coefficient "4" in Eq. d52l is the classical value, i. e. for 
each fixed even N the coefficient D(N, s) approaches 4 with increasing s ll79l . 
The approximate spectrum, (1521. is similar to that of two spins, Sa and Sb, each of spin quan- 
tum number Ns/2, that are coupled by an effective interaction of strength AJ/N. Therefore, 
one can equally well say, that the approximate rotational band considered in (l52l) is associated 
with an effective Hamilton operator 

4 J 



^-approx 



N 



S 1 



n2 



r?2 



(53) 



where the two sublattice spins, Sa, Sb, assume their maximal value Sa = S B = Ns/2. Hamil- 
tonian d53t is also known as Hamiltonian of the Lieb-Mattis model which describes a system 
where each spin of one sublattice interacts with every spin of the other sublattice with equal 
strength lT72"llg4"l . 
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Fig. 12: The low-lying levels of a spin ring, N = 6 and s = 5/2 in this example, can be grouped 
into the lowest (Lande) band, the first excited (Excitation) band and the quasi-continuum ( QC). 
For the spin levels of the L- and E-band k is given in brackets followed by the energy. Arrows in- 
dicate strong transitions from the L-band. Associated numbers give the total oscillator strength 
fofor these transitions. With friendly permission by Oliver Waldmann [81 /. 

It is worth noting that this Hamiltonian reproduces more than the lowest levels in each subspace 
H(S). At least for bipartite systems also a second band is accurately reproduced as well as 
the gap to the quasi-continuum above, compare Figure [121 and Ref. ED- This property is 
very useful since the approximate Hamiltonian allows the computation of several observables 
without diagonalizing the full Hamiltonian. 

It is of course of utmost importance whether the band structure given by the approximate Hamil- 
tonian d53t persists in the case of frustrated molecules. It seems that at least the minimal ener- 
gies still form a rotational band which is understandable at least for larger spin quantum num- 
bers s taking into account that the parabolic dependence of the minimal energies on S mainly 
reflects the classical limit for a wide class of spin systems [85|. 

The following example demonstrates that even in the case of the highly frustrated molecule 
{Mo 72 Fe3o} the minimal energies arrange as a "rotational band" 7 . In the case of {Mo 72 Fe 3 o} 

7 Work done with Matthias Exler, Universitat Osnabriick, Germany. 
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the spin system is decomposable into three sub-lattices with sub-lattice spin quantum numbers 
Sa, Sb, and Sc ll79l [80l. The corresponding approximate Hamilton operator reads 



-n approx — O 



(54) 



where S is the total spin operator and the others are sub-lattice spin operators. D and 7 are 
allowed to deviate from their respective classical values, D = 6 and 7 = 1, in order to correct 
for finite s. 




Fig. 13: DMRG eigenvalues and lowest rotational band of the s = 5/2 icosidodecahedron; 
m = 60 was used except for the lowest and first exited level which were calculated with m = 
120. 



We use the DMRG method to approximate the lowest energy eigenvalues of the full Hamiltonian 
and compare them to those predicted by the rotational band hypothesis (l54l . Fig. [13] shows the 
results and a fit to the lowest rotational band. Assuming the same dependence on m as in 
the s — 1/2 case, the relative error of the DMRG data should also be less than 1%. The 
agreement between the DMRG energy levels and the predicted quadratic dependence is very 
good. Nevertheless, it remains an open question whether higher lying bands are present in such 
a highly frustrated compound. 



Magnetization jumps 




Fig. 14: Structure of the icosidodecahedron (l.h.s.) and the cuboctahedron (r.h.s.). 

Although the spectra of many magnetic molecules possess a rotational band of minimal energies 
E m in(S) and although in the classical limit, where the single-spin quantum number s goes to 
infinity, the function E min (S) is even an exact parabola if the system has co-planar ground 
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states [85|, we 8 find that for certain coupling topologies, including the cuboctahedron and the 
icosidodecahedron (see Fig. PHI), that this rule is violated for high total spins ll67l l86l . More 
precisely, for the icosidodecahedron the last four points of the graph of E min versus S, i. e. the 
points with S = S max to S = S max — 3, lie on a straight line 

E min (S) = mJs 2 - 6Js(30s - S) . (55) 

An analogous statement holds for the last three points of the corresponding graph for the cuboc- 
tahedron. These findings are based on numerical calculations of the minimal energies for several 
s both for the icosidodecahedron as well as for the cuboctahedron. For both and other systems 
a rigorous proof of the high spin anomaly can be given 115711571. 

The idea of the proof can be summarized as follows: A necessary condition for the anomaly is 
certainly that the minimal energy in the one-magnon space is degenerate. Therefore, localized 
one-magnon states can be constructed which are also of minimal energy. When placing a second 
localized magnon on the spin array there will be a chance that it does not interact with the 
first one if a large enough separation can be achieved. This new two-magnon state is likely 
the state of minimal energy in the two-magnon Hilbert space because for antiferromagnetic 
interaction two-magnon bound states do not exist. This procedure can be continued until no 
further independent magnon can be placed on the spin array. In a sense the system behaves as if 
it consists of non-interacting bosons which, up to a limiting number, can condense into a single- 
particle ground state. In more mathematical terms: In order to prove the high-spin anomaly 
one first shows an inequality which says that all points (S, E min (S)) lie above or on the line 
connecting the last two points. For specific systems as those mentioned above what remains to 
be done is to construct particular states which exactly assume the values of E min corresponding 
to the points lying on the bounding line, then these states are automatically states of minimal 
energy. 




2 4 6 8 10 12 U 5 1 15 2 2 5 3.0 3.5 



S g MB B/J 

Fig. 15: Icosidodecahedron: L.h.s. - minimal energy levels E min (S) as a function of total spin 
S. R.h.s. - magnetization curve atT = H67V . 

The observed anomaly - linear instead of parabolic dependence - results in a corresponding 
jump of the magnetization curve M. versus B, see Fig. [15] In contrast, for systems which obey 
the Lande interval rule the magnetization curve at very low temperatures is a staircase with 
equal steps up to the highest magnetization. The anomaly could indeed be observed in magneti- 
zation measurements of the Keplerate molecules {Mo72Fe3o}. Unfortunately, the magnetization 
measurements ll36ll80l performed so far suffer from too high temperatures which smear out the 
anomaly. 

8 Work done with Heinz-Jiirgen Schmidt, Universitat Osnabriick, Andreas Honecker, Universitat Braunschweig, 
Johannes Richter and Jorg Schulenburg, Universitat Magdeburg, Germany. 
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Nevertheless, it may be possible to observe truly giant magnetization jumps in certain two- 
dimensional spin systems which possess a suitable coupling (e. g. Kagome) [86|. In such 
systems the magnetization jump can be of the same order as the number of spins, i. e. the jump 
remains finite - or in other words is macroscopic - in the thermodynamic limit N — ► oo. Thus, 
this effect is a true macroscopic quantum effect. 

4 Dynamics 

In this section I would like to outline two branches - tunneling and relaxation - where the 
dynamics of magnetic molecules is investigated. The section is kept rather introductory since 
the field is rapidly evolving and it is too early to draw a final picture on all the details of the 
involved processes. 

4.1 Tunneling 

Tunneling dynamics has been one of the corner stones in molecular magnetism since its very 
early days, see e. g. (EHUD M M H3- 

The subject can roughly be divided into two parts, one deals with tunneling processes of the 
magnetization in molecules possessing a high ground state spin and an anisotropy barrier, the 
second deals with the remaining tunneling processes, e. g. in molecules which have an S = 
ground state. 




M M 



Fig. 16: Sketch of the tunneling barrier for a high spin molecule with S = 10, l.h.s. without 
magnetic field, r.h.s. with magnetic field, compare Eq. 0. The arrows indicate a possible 
resonant tunneling process. 

As already mentioned in section I3TTI some molecules like Mn 12 and Fe 8 possess a high ground 
state spin. Since the higher lying levels are well separated from the low-lying S = 10 levels 
a single-spin Hamiltonian ©, which includes an anisotropy term, is appropriate. Figure \TG\ 
sketches the energy landscape for an anisotropy term which is quadratic in S z . If the Hamilto- 
nian includes terms like a magnetic field in x-direction that do not commute with S z resonant 
tunneling is observed between states | S, M ) and | S, —M). This behavior is depicted on the 
l.h.s. of Fig.[J3for the transition between M = —10 and M = 10. If an additional magnetic 
field is applied in ^-direction the quadratic barrier acquires an additional linear Zeeman term 
and is changed like depicted on the r.h.s. of Fig. [T5J Now tunneling is possible between states 
of different \M\, see e. g. EE). 

It is rather simple to model the tunneling process in the model Hilbert space of S = 10. i. e. 
a space with dimension 2S + 1 = 21. Nevertheless, in a real substance the tunneling process 
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is accompanied and modified by other influences. The first major factor is temperature which 
may enhance the process, this leads to thermally assisted tunneling [92|. Each such substance 
hosts phonons which modify the tunneling process, too, resulting in phonon assisted tunneling 
11951 1931 1951 l9"5l . Then local dipolar fields and nuclear hyperfine fields may strongly affect the 
relaxation in the tunneling regime 0. In addition there may be topological quenching due to 
the symmetry of the material Il97ll98ll99ll . And last but not least describing such complicated 
molecules not in effective single-spin models but in many-spin models is still in an unsatisfac- 
tory state, compare [100]. 




Fig. 17: Sketch of the tunneling process between N eel-like states on a spin ring. Without loss of 
generality the state on the l.h.s. will be denoted by \ Neel, 1 ) and the state on the r.h.s. will be 
denoted by \ Neel, 2 ). 

Another kind of tunneling is considered for Heisenberg spin rings with uniaxial single-ion 
anisotropy. Classically the ground state of even rings like Na:Fe 6 and Cs:Fes is given by a 
sequence of spin up and down like in Fig.[T7l It now turns out that such a Neel-like state, which 
is formulated in terms of spin-coherent states d40l . contributes dominantly to the true ground 
state as well as to the first excited state if the anisotropy is large enough [69|. Thus it is found 
that the ground state | E ) and the first excited state | Ei ) can be approximated as 

\E ) « -=(|N6el,l)± | Neel, 2)) (56) 
V2 

\Ei) « i=(|Neel,l) T | Neel, 2 » , 

where the upper sign is appropriate for rings where the number of spins N is a multiple of 4, 
e. g. N = 8, and the lower sign is for all other even N. 

Therefore, the tunneling frequency is approximately given by the gap between ground and first 
excited state. Experimentally, such a tunnel process is hard to observe, especially since ESR is 
sensitive only to the total spin. What would be needed is a local probe like NMR. This could be 
accomplished by replacing one of the iron ions by another isotope. 

The tunneling process was further analyzed for various values of the uniaxial single-ion anisotropy 
1 101 1. Since in such a case the cyclic shift symmetry persists, k is still a good quantum number. 
Therefore, mixing of states is only allowed between states with the same k quantum number. 
This leads to the conclusion that the low-temperature tunneling phenomena can be understood 
as the tunneling of the spin vector between different rotational modes with AS = 2, compare 
Fig. Eland the subsection on rotational bands on pagefT9l 

4.2 Relaxation dynamics 

In a time-dependent magnetic field the magnetization tries to follow the field. Looking at this 
process from a microscopic point of view, one realizes that, if the Hamiltonian would commute 
with the Zeeman term, no transitions would occur, and the magnetization would not change a 
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Fig. 18: Energy spectrum of spin rings with N = 6 and vanishing anisotropy at two magnetic 
fields drawn as a function of the magnetic quantum number M. The dashed curves represent the 
lowest-lying parabolas E min (M) discussed in section \TM A white or black circle indicates that 
a state belongs to k = or k = N/2, compare Fig. \T7\ States belonging to one spin multiplet 
are located on straight lines like that plotted in panel (a) for the 5 = 4 multiplet. With friendly 
permissions by Oliver Waldmann 1 101 J. 



tiny bit. There are basically two sources which permit transitions: non-commuting parts in the 
spin Hamiltonian and interactions with the surrounding. In the latter case the interaction with 
phonons seems to be most important. 

Since a complete diagonalization of the full Hamiltonian including non-commuting terms as 
well as interactions like the spin-phonon interaction is practically impossible, both phenomena 
are modeled with the help of rate equations. 



B 



Fig. 19: Schematic energy spectrum in the vicinity of an avoided level crossing. The formula by 
Landau, Zener, and Stiickelberg (1571 ) approximates the probability p for the tunneling process 
from Mi to M 2 . 

If we start with a Hamiltonian H , which may be the Heisenberg Hamiltonian and add a non- 
commuting term H' like anisotropy the eigenstates of the full Hamiltonian are superpositions of 
those of H . This expresses itself in avoided level crossings where the spectrum of H would 
show level crossings, compare Fig.[l3 Transitions between eigenstates of H , which may have 
good M quantum numbers, can then effectively be modeled with the formula by Landau, Zener, 
and Stiickelberg lfTMfT0^fTMfT05lfT06l 

VTA 2 

p = 1 - exp * ■ ; — -j — } . (57) 

• ^ 2hg f i B \M 1 -M 2 \ f t B ( 



A denotes the energy gap at the avoided level crossing. 
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The effect of phonons is taken into account by means of two principles: detailed balance, which 
models the desire of the system to reach thermal equilibrium and energy conservation, which 
takes into account that the energy released or absorbed by the spin system must be absorbed 
or released by the phonon system and finally exchanged with the thermostat. The interesting 
effects arise since the number of phonons is very limited at temperatures in the Kelvin-range 
or below, thus they may easily be used up after a short time (phonon bottleneck) and have to 
be provided by the thermostat around which needs a characteristic relaxation time. Since this 
all happens in a time-dependent magnetic field, the Zeeman splittings change all the time and 
phonons of different frequency are involved at each time step. In addition the temperature of 
the spin system changes during the process because the equilibration with the thermostat (liq. 
Helium) is not instantaneous. More accurately the process is not in equilibrium at all, especially 
for multi-level spin systems. Only for two-level systems the time-dependent occupation can be 
translated into an apparent temperature. In essence the retarded dynamics leads to distinct hys- 
teresis loops which have the shape of a butterfly H107Lll08lll09l . For more detailed information 
on dissipative two-level systems the interested reader is referred to Ref. HI 10111091 . 



5 Magnetocalorics 

The mean (internal) energy, the magnetization and the magnetic field are thermodynamic ob- 
servables just as pressure and volume. Therefore, we can design thermodynamic processes 
which work with magnetic materials as a medium. This has of course already been done for 
a long time. The most prominent application is magnetization cooling which is mainly used 
to reach sub-Kelvin temperatures 11111 . The first observation of sub-Kelvin temperatures is a 
nice example of how short an article can be to win the Nobel prize (Giauque, Chemistry, 1949). 
Meanwhile magnetization cooling is used in ordinary refrigerators. 



768 LETTERS TO THE EDITOR 

Attainment of Temperatures Below 1° Absolute T>y Demagnetization of GdjfSOJj.SHjO 



We have recently carried out some preliminary experi- 
ments on the adiabatic demagnetization of Gds(SO.|)s 

■ 8H;G at the temperatures of liquid helium. As previously 
predicted by one of us r a large fractional lowering of the 
absolute temperature was obtained. 

An iron-free solenoid producing a field of about 8000 
gauss was used for all the measurements. The amount of 
Gd^SG-sXvSltaQ was 61 g. The observations were checked 
by many repetitions of the cooling. The temperatures were 
measured by means of the inductance of a coil surrounding 
the gadolinium sulfide. The coil was immersed in liquid 
helium and isolated from the gadolinium by means of an 
evacuated space. The thermometer was in excellent 
agreement with the temperature of liquid helium as 
indicated by its vapor pressure down to 1.5°K. 



On March 19 3 starting at a temperature of about 3.4°K, 
the material cooled to G.53K. On April 8, starting at 
about 2", a temperature of 0..M°K was reached. On 
April 9, starting at about 1.5°, a temperature of 0.25°K 
was attained. 

Ft is apparent that it will be possible to obtain much 
lower temperatures, especially -when successive demagneti- 
zations are utilized. 

W. F. GlAUQUE 

D. P. MacDougall 

Department gf Chemistry, 
University of California, 
Berkeley, California, 
April 12, 1933. 



Fig. 20: The first observation of sub-Kelvin temperatures Ml IV is a nice example of how short 
an article can be to win the Nobel prize (Giauque, Chemistry, 1949). 



In early magnetocaloric experiments simple refrigerants like paramagnetic salts have been used. 
We will therefore consider such examples first. 

For a paramagnet the Hamiltonian consists of the Zeeman term only. We then obtain for the 



X5.28 



Jiirgen Schnack 



partition function 



Then the magnetization is 

M(T, B, N) = Ngfi B {{s + l/2)coth{(3gfi B B(s + 1/2)] - l/2sinh[/%i BJ B/2]} , (59) 
and the entropy reads 

- M 'lBBH -^^- (60) 

Besides their statistical definition both quantities follow from the general thermodynamic rela- 
tionship 

( d F\ ( d F\ 

dF = — dT+ — dB = -SdT - MdB , (61) 



y dT J B \dB /T 
where F(T, B, N) = -k B T ln[Z(T, B, N)]. 

Looking at Eq. (I5HI it is obvious that all thermodynamic observables for a paramagnet depend 
on temperature and field via the combination B/T, and so does the entropy. Therefore, an 
adiabatic demagnetization (S = const) means that the ratio B/T has to remain constant, and 
thus temperature shrinks linearly with field, i.e. 

d T \ ^ £tra T 

(62) 



dB) s B 

This situation changes completely for an interacting spin system. Depending on the interactions 
the adiabatic cooling rate |^ can be smaller or bigger than the paramagnetic one (l62l) and even 
change sign, i.e. one would observe heating during demagnetization. It is nowadays understood 
that the cooling rate acquires extreme values close to phase transitions due to the excess entropy 
associated with such processes [ 1121I1131I114U115I . 




Fig. 21: L.h.s.: Sketch of an antiferromagnetically coupled spin dimer. R.h.s.: Dependence of 
the energy levels on magnetic field for an antiferromagnetically coupled spin-1/2 dimer. At the 
critical field B c the lowest triplet level crosses the ground state level with S = 0. 



In the following this statement will be made clear by discussing the example of a simple anti- 
ferromagnetically coupled spin-1/2 dimer (Tig. I2T1 l.h.s.). In a magnetic field such a system 
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Fig. 22: L.h.s.: Entropy of the dimer (Fis. \21l as function of B and T. R.h.s.: Heat capacity of 
the dimer as function of B and T. 



experiences a "quantum phase transition" if the lowest triplet level crosses the original ground 
state with S = 0, see Fig. EH r.h.s.. Although one would hesitate to call such an ordinary ground 
state level crossing quantum phase transition it nevertheless is one. At T = the magnetization 
Ai(T = 0, B) is a non-analytic function of the magnetic field B. At the critical field B c where 
the levels cross the magnetization exhibits a step. 

In addition the entropy, which at T = is zero for the non-degenerate ground state acquires a 
finite value at the critical field B c due to the degeneracy of the crossing levels. This enhancement 
remains present even at temperatures T > 0, compare Fig. E21 l.h.s.. In addition the heat 
capacity varies strongly around the critical field as is shown in Fig.|22l r.h.s.. 



T-B isentropes 




0.5 1 1.5 2 2.5 3 



Fig. 23: Isentrops of the spin dimer. The straight lines show the behavior of a paramagnet for 
comparison. B is along the x-axis, T along the y-axis. 



The behavior of the entropy as well as of the heat capacity explains how the adiabatic cooling 
rate 

dT\ t {§Mt (63) 



dBj s C 

depends on field and temperature. Figure |23] shows the isentropes of the antiferromagnetically 
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coupled dimer both as function of field B and temperature T. The straight lines show the 
behavior of a paramagnet for comparison. Three regions are highlighted. 

• a: For low fields and high temperatures Sjs is smaller than for the paramagnet. 



• b: For high fields and high temperatures the interacting system assumes the paramagnetic 
limit, i.e. |^ is the same in both systems. 

• c: For low temperatures and fields just above the critical field is much bigger than the 
cooling rate of the paramagnet. 

• Not highlighted but nevertheless very exciting is the region at low temperature just below 
the critical field where the "cooling" rate |^ has the opposite sign, i.e. upon demagnetiz- 
ing the system heats up and upon magnetizing the system cools down. 

The rate |^ d63l depends directly on the derivative of the entropy with respect to the magnetic 
field. Therefore, it is clear that the effect will be enhanced if a high degeneracy can be ob- 
tained at some critical field. This is indeed possible in several frustrated materials where giant 
magnetization jumps to saturation are observed I86l lll6llll41lll5l . 



cuboctahedron, s=1/2 



ring, N=12, s=1/2 





0.0 0.2 O.i 6 8 1.0 1.2 1.4 16 

B/B sat 



Fig. 24: Isentropes of an antiferromagnetically coupled spin-1/2 cuboctahedron (l.h.s.) in 
comparison to the isentropes of an antiferromagnetically coupled spin ring of 12 spins with 
s = 1/2 (r.h.s.). The structure of the cuboctahedron is shown in Fis.\14\(r.h.s. ). 



The following Figure l24l compares isentropes of an antiferromagnetically coupled spin-1/2 
cuboctahedron (l.h.s.) with the isentropes of an antiferromagnetically coupled spin ring of 12 
spins with s — 1/2 (r.h.s.). It is clearly visible that the cuboctahedral spin system has steeper 
isentropes than the ring system, and the reason is that the cuboctahedron features a larger jump 
to magnetization saturation (AM = 2) than the spin ring (AM = 1) which is connected with 
an enhanced degeneracy and thus higher entropy. 

Summarizing, one can say that low-dimensional frustrated spin systems and especially magnetic 
molecules are substances with an interesting magnetocaloric behavior and may turn out to be 
useful new refrigerants for special applications. 
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